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Abstract 

We  regard  the  transition  layer  near  the  bound  of  inhomogeneous  media  on  the  base  of  the 
introduced  here  concept  of  the  generalized  wave  modes  (GWM).  We  show  in  present  work  that 
previously  described  one-dimensional  nonlocal  model  of  inhomogeneous  media  [1]  is  a two- wave 
approximation  for  electrodynamic  properties  of  the  medium.  The  second,  evanescent,  wave  in  the 
model  represents  a sum  of  wide  spectrum  of  GWM  and  can  be  used  for  approximate  description  of 
the  transition  layer. 


1.  Introduction 

The  transition  layer  is  relatively  thin  area  near  the  bound  of  medium,  where  electrodynamic  property 
cannot  be  described  adequately  in  terms  of  £,fl  parameters.  Usually,  inhomogeneous  media,  such  as 
composite  materials,  have  a depth  of  the  transition  layer  approximately  the  same  as  a characteristic 
size  a of  the  structure.  But  one  can  expect  that  in  special  cases  this  depth  can  be  considerably  larger 
than  a,  at  least  for  the  materials  with  sufficiently  strong  interaction  between  the  inclusions. 

Taking  into  account  the  presence  of  the  transition  layer  for  inhomogeneous  materials,  one  can 

provide  the  general  conclusion  that  the  measured  effective  parameters  £e , fle  of  them  must  always 
differ  from  those  of  £e,  fle,  which  one  can  calculate  for  the  infinite  (boundless)  medium. 

The  most  transparent  example  for  this  conclusion  can  be  the  finely  stratified  medium  (the  Rytov 
solution  [2]),  which  contain  dielectric  and  metal  layers.  Let  us  consider  only  TM  modes  in  it.  There 
exist  only  one  propagating  (principal)  wave  mode  in  such  a structure  and  all  higher  order  modes  are 
evanescent  in  the  quasi-static  case.  Acher  et  al  [3]  used  formulas  for  the  £e,  fle  parameters,  which 
correspond  to  the  principal  wave  in  such  a structure,  for  calculation  the  intrinsic  permeability  of  thin 
metal  films.  In  the  limit  of  strong  skin  effect  those  formulas  give  the  result  £e  =£,  He  = 1 , which  is 

not  correct,  because  do  not  reflects  the  diamagnetic  properties  of  such  a structure.  Our  theoretical 
analysis  and  the  experimental  examination  permit  us  to  found  the  correct  formulas  for  the  considered 
case: 


e l e d 

£ =£->£,  A*  =T<1  (a) 

a l 

where  l is  a period  of  the  layers  and  d is  the  width  of  dielectric  layers.  The  difference  between  the  first 
and  the  second  formulas  can  be  explained  by  the  fact  of  excitation  the  wide  spectrum  of  evanescent 
modes  near  the  plane  boundary  of  this  medium.  These  modes  form  the  transition  layer,  which  depth  is 
close  to  the  skin  depth.  Therefore,  the  formulas  in  [3]  are  not  correct,  if  the  skin  effect  is  noticeable. 
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2.  Determination  of  the  Generalized  Wave  Modes  in  Inhomogeneous 
Medium 

A 

Let  us  consider  arbitrary  inhomogeneous,  but  statistically  homogeneous  and  isotropic  in  x and  y 
directions  medium.  It  is  conveniently  to  regard  the  transition  layer  on  the  boundary  on  the  base  of 
conception  of  the  generalized  wave  modes  (GWM).  Such  a mode  we  determine  as  a generalized  wave  \ 

process  in  the  medium,  when  the  spatially  mean  constituents  Eo  and  Ho  of  the  electromagnetic  field  in 
it  spread  as  a plane  wave  (factor  exp (icot)  is  omitted) 

E0(z)  = £060exp(-ifcz),  H0(z)  = ff0fi0exp(-/*z),  fi0=[2xfi0]  (1) 

This  definition  is  consistent  with  the  theory  of  the  space  dispersion  [4]  where  only  spatially  means 
constituents  of  the  field  are  taken  into  consideration.  In  present  work  we  shall  take  into  account  the 
spatially  inhomogeneous  constituents  of  the  field  also.  For  this  purpose  we  provide  decomposition  of 
tangential  components  of  the  whole  electromagnetic  field  in  arbitrary  plane  (x,  y)  in  form  of  finite 
series  of  the  spatial  harmonics  (SH).  For  example,  it  can  be  done  with  help  of  2D  FFT,  but  we  suppose 
that  each  SH  is  a sum  of  both  elementary  plane  waves  in  the  (x,  y)  plane  with  the  same  wavenumber. 

Further  we  shall  apply  one-dimension,  rather  then  two-dimensions  indexing  of  SH  for  the  simplicity. 

Providing  transformations  of  the  Maxwell's  equations  by  analogy  to  the  method,  have  been  described 
in  [1],  one  can  obtain  the  followed  system  of  equations 


N-l 

kEn  = (Qfi  o nmH „ 

m- 0 
AM 

kHn  =coe0  Yj£  nmEm 

m- 0 


n — 0,...,N -l 


(2) 


Here  we  use  SI  system,  Ea,  Hn  are  the  amplitudes  of  the  spatial  harmonics,  N - the  number  of  the 
harmonics,  enm,  flnm  are  complex  parameters,  which  depend  on  £,  fl  distribution  in  the  (x,  y)  plane 

and  polarization  of  the  SH  only.  Naturally,  the  zero  index  corresponds  to  the  spatially  homogeneous 
constituents  of  the  field.  One  can  show  with  help  of  the  Maxwell's  equations  that  for  the  defined 
polarization  of  the  E0,  Ho  constituents,  polarization  of  the  fields  of  each  SH  in  the  defined  (x,  y)  plane 
do  not  depend  on  the  wavenumber  k of  the  mode  (1).  Hence  the  same  is  true  for  the  parameters 
£nm , Hnm . Moreover,  it  seems  natural  that  those  parameters  do  not  depend  on  z for  statistically 

homogeneous  media,  although  we  have  no  proof  of  it  at  present  time.  Nevertheless,  if  there  are  some 
fluctuations  of  those  parameters  in  z direction,  we  always  can  take  the  mean  values  of  them. 

Therefore,  the  system  (2)  is  a linear  system  of  equations.  It  can  be  rewritten  in  the  equivalent  and 
a more  convenient  form,  as  a combination  of  two  matrix  equations: 


N- 1 


A*E  = n2-E,  A = [a„ J,  anm  = JX*  •£*,,  E = (£0,£„.. 


*=0 

AM 


B*H=n2  -Hj  £ = [&„„1  -Ab.,  H = 


(3) 


k= 0 


Here  n = k/k0  is  a refractive  index  of  the  GWM,  which  can  be  obtained  as  a solution  to  the  two 
equivalent  dispersion  equation 


'tw 


A 


det(^  - n 2 /)=  0 det(g  - « 2 /)=  0 


(4) 


413 


where  / is  a unite  matrix.  It  is  clearly  that  the  eigenvectors  E,H_oi  the  corresponding  matrixes  A,  B 

describe  the  spatially  spectrum  of  SH  for  respective  GWM.  Therefore,  in  contradistinction  to  the 
conventional  concept  of  wave  mode,  the  GWM,  travelling  trough  the  inhomogeneous  medium,  keeps 
only  the  spatial  spectrum  of  transversal  distribution  of  the  electromagnetic  field. 

If  we  pass  to  the  limit  of  the  continuous  spatially  spectrum,  the  characteristic  polynomial 

det^  ~n2l)  transfer  to  an  entire  function  [5],  which  have  isolated  zeros  only.  Hence  the  spectrum  of 

GWM  for  the  medium,  which  is  described  by  the  equation  (4),  is  always  discreet.  It  is  clear  that 
similar  to  the  finely  stratified  medium  [2],  the  only  principal  GWM  is  usually  propagating  and  all 
higher  order  modes  are  evanescent  in  quasi-static  case.  Now  we  can  conclude  that  by  analogy  to  the 
Rytov  solution,  the  transition  layer  is  formed  by  a series  of  evanescent  generalized  modes,  which  are 
excited  by  an  incident  wave. 


3.  Two- Wave  Approach  to  Description  of  the  Transition  Layer 

Let  us  show  that  one-dimensional  nonlocal  model  of  the  inhomogeneous  medium  [1]  corresponds  to  a 
two-wave  approximation  for  electrodynamic  description  of  the  medium.  For  the  uniformity,  we 
rename  the  8 constitutive  parameters  introduced  in  [1]  to  e^,  Hkm  • Then  we  can  present  4 equations 
obtained  in  [1]  in  the  next  form 

k ■ E0  =cofi()(p.Q0H0+  k ■ H0  = (Q£0(£q0E0  + eoxEx) 

k-E)  =(Qfi0(jll0H0  + fluHx)  k-Hx  = CQ£0(£X0E0+£XXEX) 

Here  E\  and  Hx  are  the  effective  complex  amplitudes  of  the  spatially  inhomogeneous  constituents  of 
the  tangential  components  of  the  field  [1].  These  amplitudes  can  be  express  in  the  amplitudes  of 
spatial  harmonics  (2) 

I AM  „ /AM  “ 

*,=J lEl,  (6) 

V m-\  V m= 1 

Parameters  £km , jlkm  can  be  calculated  by  formulas,  which  followed  from  [1],  where  functions  of 

distribution  of  the  inhomogeneous  constituents  of  the  field  correspond  to  the  principal  wave  in  the 
media. 

We  can  rewrite  the  system  (5)  in  form  of  two  systems  by  analogy  to  eq.(3): 

\aoo^o+ao\^i =n  -Eo  lb00H0+b0lHx=n  -H0  ^ 

laio-^o  ^a\\E\  — n ' E\  **” b\\Hx  = n 'Hx 

Here  coefficients  aksn  and  felon  are  depended  on  parameters  £km , in  similar  way  to  (3): 

aoo  =^ooeoo  + j“oi£,o> 
aio  =iUi0£oo+Mi,£io> 

^oo  = eoo  ^oo  + eoi  /*io  > 

^10  =eio/ioo  + £11^10> 

Solution  of  Eq.(7)  can  be  obtained  by  analogy  to  Eq.  (4)  through  dispersion  equations 


aoi  =^ooeoi  + A*oien 
a\\  =Ml0e01  + ^llell 
^01  =e00^01  +£01^11 
^11  =eioMoi  +£11^11 
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detQ2-«2/2)=0<=>det^2-n2/2)=0  (9) 

where  index  2 corresponds  to  2x2  matrixes  of  Eq.(7).  Hence  the  equations  (9)  have  two  solutions  for 
the  value  n.  Naturally,  the  first  root  of  (9)  must  correspond  to  the  principal  wave  in  the  medium.  The 
second  solution  represents  an  effective  GMW,  which  is  a superposition  of  all  higher  evanescent 
generalized  modes  in  it.  Therefore  one  can  expect  that  the  second  wave  of  the  one-dimensional  model 
can  be  used  for  the  approximate  description  of  the  transition  layer  of  inhomogeneous  media.  We  mean 
here  the  case  of  plane  boundary  of  the  medium  and  plane  wave,  which  normally  fall  on  it. 

The  model  of  medium,  where  two  wave  modes  can  be  excited,  are  well-known  from  theory  of  the 
spatial  dispersion  [4],  although  only  propagating  waves  was  taken  into  consideration  there.  The 
principal  point  of  this  problem  is  to  find  the  additional  bound  condition,  which  is  necessary  for  sewing 
together  tangential  components  of  electromagnetic  fields  on  the  boundary  for  incident  wave  and  both 
modes.  The  same  problem  arises,  when  we  consider  the  transition  layer  in  two-wave  approximation. 
As  an  additional  boundary  condition,  we  can  use  in  this  case  the  continuity  of  tangential  components 
of  the  spatial  inhomogeneous  constituent  on  the  plane  boundary.  This  requirement  is  equivalent  to  the 
introducing  of  the  complex  effective  characteristic  impedance  Zoi  for  the  spatial  inhomogeneous 
constituent  of  the  field,  which  spread  outside  the  bound.  Then,  if  the  bound  lies  at  z = 0,  we  can  write 

the  additional  bound  condition  in  the  form  of  El  (0)/H{  (0)  = Z01 . The  value  of  Z01  = Zq-J 
can  be  obtained  from  the  formulas  in  [1],  where  one  must  use  values  £ = 1,  jU  = 1 for  the  free  space. 

We  have  examined  the  two-wave  approach  in  particular  for  the  case  of  the  finely  stratified  medium 
and  get  good  results.  For  example,  if  such  a structure  contains  both  dielectric  and  metal  layers,  we 
have  the  correct  result  (a)  for  the  limit  of  strong  skin  effect  in  metal  layers. 


4.  Conclusion 

The  transition  layer  in  inhomogeneous  media  can  be  described  with  help  of  series  of  the  GWM.  But 
we  think  that  the  simple  two-wave  approach  for  this  problem  is  quite  adequate  in  many  cases.  Such  an 
approach  is  seemed  more  acceptable,  especially  in  case,  when  we  investigate  both  theoretically  and 
experimentally  the  influence  of  the  transition  layer  on  electromagnetic  properties  of  complex 
composite  materials. 
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